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ABSTRACT 

Context. One of the most dramatic manifestations of solar activity are large-scale coronal bright fronts (CBFs) observed in extreme 
ultraviolet (EUV) images of the solar atmosphere. To date, the energetics and kinematics of CBFs remain poorly understood, due to 
the low image cadence and sensitivity of previous EUV imagers and the limited methods used to extract the features. 
Aims. In this paper, the trajectory and morphology of CBFs was determined in order to investigate the varying properties of a sample 
of CBFs, including their kinematics and pulse shape, dispersion, and dissipation. 

Methods. We have developed a semi-automatic intensity profiling technique to extract the morphology and accurate positions of 
CBFs in 2.5-10 min cadence images from STEREOfEVYl. The technique was applied to sequences of 171 A and 195 A images from 
STEREO/EUVl in order to measure the wave properties of four separate CBF events. 

Results. Following launch at velocities of ~240-450 km s -1 each of the four events studied showed significant negative acceleration 

ranging from 290 to -60 m s~ 2 . The CBF spatial and temporal widths were found to increase from ~50 Mm to ~200 Mm and 

~100 s to ~1500 s respectively, suggesting that they are dispersive in nature. The variation in position-angle averaged pulse-integrated 
intensity with propagation shows no clear trend across the four events studied. These results are most consistent with CBFs being 
dispersive magnetoacoustic waves. 
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1. Introduction 

"EIT waves" were first observed in the solar corona using the 
SOlar and Heliospheric Observ atory (SOi/0)/Ext reme ultra- 
violet Imaging Te lescope (EIT; Moses et"afl Il997l) and anal- 
ysed in detail by iThompson et alJ dl998l) . Thev have been a 
source of much controversy and debate in the solar physics 
community since this initial observation, with different authors 
sugge sting that the y are alternatively fast-mode MHD waves 



(e.g.: IWangj 



IVeronig et al. 



2000|: IWarmufh et all l2004at Long etail l2008t 
iGope " 



2008 



iGopalswamv et all l2009l) . a result of the re- 



structuring of the ma gnetic field during the eruption of a coronal 
mass ejection (CME; [Chen et al.ll2002l l2005b lAttrill et al.ll2006l 
l2007t iDelannee et alJl2007l l2008h or a coronal MHD soliton 
dWiils -Davev et al.ll2007t) . It should be noted at this point that 
the name "EIT wave" is typically used for historical reasons. 
To reflect the uncertainty surrounding the physical interpreta- 
tio n of "EIT waves" we shal l adopt the suggested nomenclature 
of iGallagher & Long! (EoToh and refer to these disturbances as 
coronal bright fronts (CBFs) throughout this paper. 

The uncertainty in the nature of CBFs has arisen from con- 
flicting results being drawn from the same observations. A 
pseudo-wave interpretation was proposed following observa - 
tions of stationary bright fronts (Delanne e & Aulanierlll999l) . 
a strong correlation with CMEs dBiesecker et al.l |2002[) and a 
lower than expected estimated pulse velocity dwills-Davev et alJ 
20Q7h. However, observa t ions of refraction (IWangl I20001 



Ofm an & Thompsonl l2002t IVeronig et alJ I2006I) and reflection 
(IGopalswamv et aD I2009I) of CBFs at coronal hole boundaries 
would appear to suggest a wave interpretation. A conclusive re- 
sult has been hampered by the diffuse nature of CBFs and the 



relatively low temporal cadence of the observing instruments, 
both of which make it difficult to characterize their true nature. A 
full review of CBFs, their morphology, kinematics, relationship 
to other so lar phenomena, and theoreti cal i nterpretations may 
be fou nd in lWills-Davev & Attrilll (l2010h and lGallagher & Longl 
J20TOI) . 

The first observations of CBFs were made using EIT 
with an effective cadence of ~12 minutes in the 195 A 
passband. Initial estimates of the kinematics of these dis- 
turbances, using a point-and-click methodology applied to 
running-differ ence images, estimated th e average velocity at 
~ J 89 km s _1 (IThompson & Mversll2009l). A higher velo city of 
311 ± 111 km s" 1 was found by IWarmufh et all (l2004al) using 
additional passbands to compensate for the lack of 195 A im- 
ages. Thi s is comparable to the ra nge of Alfven speeds pre- 
dicted by IWills-Davev et alJ (120071 -215-1500 k m s" 1 ). The 
Solar TErrestrial RElations Observatory (SrE7?E6>: lKaiser et al.l 
120081) mission with the Extreme Ultra Violet Imager (EUVI; 
IWuelser et al.l20 04) instrument has led to new results. EUVI has 
an effective observing cadence of up to 1.5 minutes (ten times 
that of EIT), allowing for an improved estimate of the kinemat- 
ics of these dis turbances. A num erical differencing technique 
was applied by lLong et all (|2P08i) to running-difference EUVI 
images to estimate a peak velocity range for a CBF of ~153 to 
475 km s _1 with the acceleration of the disturbances estimated to 
be between -413 and 816ms 2 , depending on the cad ence of the 
observations. For the same event. fVeronig et al1(l2008l) estimated 
the CBF velocity to be 460 km s _1 with an associated decelera- 
tion of -160 m s~ 2 by fitting a quadratic model t o distance-time 
measurements. The results of lLong et al.l (120081) also indicated 
that the lower cadence of SOHO/EYT had resulted in the kine- 
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matics of CB Fs being previously u nder estimated , an ob servation 
confirmed by I Veronig etafl (120081) and lMa et al.l(l2009l) . 

The kinematics of CBFs provide an insight into the true 
nature of the disturbances. Another useful physical indicator 
is the presence of pul se broadening, previously described by 
IWarmuth et all (l2004bl) . Several authors have examined CBF 
pulse broadening using observations of multiple events from dif- 
ferent pass bands, including a combination of EUV and Ha ob- 
serva tions (IWarmuth et al.ll2001t IWarmuthll20lol IVeronig et all 
1201(H) as well as observations of a disturbance in soft X- 
ray data from GOE S/SXl (IWarmuth et al.l 120051) . In contrast, 
IWills-Davevi d2006l) observed no measurable increase in the 
FWHM of a pulse using high-cadence observations of a single 
CBF across the limited field-of-view of the Transition Region 
And C oronal Explorer (TRACE). This led IWills-Davev et alj 
d2007l) to propose that CBFs were soliton-like waves which ex- 
hibit no significant dispersion with propagation. 

Observations of CBFs typically show a decreasing front in- 
tensity with p ropagation. Thi s has been no ted by many au- 
thors i ncluding IWarmuth et al.1 (120011 l2004bb and lVeronig et alj 
(l2OT0h . and is interpreted as evidence of energy conservation as 
a wave expands into a larger area. The specific geometry of the 
wave expansion can be derived by measuring the decay in pulse 
intensity and growth of pulse width, although previous efforts to 
measure these parameters accurately have been hampered by the 
small number of data p oints associated with each observed event 
( IWarmuth et al.ll2004bl) . 

In this paper we determine the kinematics of several CBFs 
and examine the variation in their width and intensity with in- 
creasing time (and hence distance). The data are presented and 
discussed in Section [2] with Section [3] detailing the analysis 
method. The results are discussed in Section |4] and some con- 
clusions drawn in Section[5] along with some thoughts on future 
work. 



2. Observations 

The data discussed here were obtained using EUVI, part of the 
Sun Earth Connection Coronal and Heliospheric Investigation 
(SECCHI; iHoward et al.l [2008b suite of instruments onboard 
the STEREO-A and STEREO-B spacecraft. EUVI is a normal- 
incidence telescope of Ritchey-Chretien design with a pixel 
scale of 1.6". It observes the Sun in four passbands (304 A, 
171 A, 195 A, and 284 A), with the peak temperature sensitiv- 
ities for each passband being approximately 0.07 MK (304 A), 
1 MK (171 A), 1.5 MK (195 A), and 2.25 MK (284 A). The 
imaging cadence of the EUVI instrument ranges from 1 .5 min- 
utes (in the 171 A passband) to 20 minutes (in the 284 A pass- 
band). 

Although EUVI takes observations in all four of these pass- 
bands, only the 171 A and 195 A passbands were used here. 
This is due to the high temporal cadence of both passbands (1.5- 
2.5 minutes for 171 A and 5-10 minutes for 195 A) and also as 
CBFs are more readily observed (i.e., of higher contrast) in these 
two passbands. CBFs have been observed in the 284 A passband 
dZhukov & Aucherdl2004l) and the 304 A passband dLong et al.l 
120081) . but the nature of the data make it difficult to use these 
passbands for more rigorous analysis. 

The CBFs were studied using de-rotated base-difference 
(BD) and percentage base-difference (PBD) images. This in- 
volved de-rotating all images for a given event to the same pre- 
event time, in order to correct for solar rotation between images, 
and then subtracting a pre-event image to produce a BD image. 
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Fig. 1. Percentage base difference images of the 2007 May 19 
(top left), 2007 Dec 07 (top right), 2009 Feb 12 (bottom left) 
and 2009 Feb 13 (bottom right) CBFs, each in the 195 A pass- 
band as seen by STEREO-B. The solid lines indicate the region 
in which the disturbance was identified using the intensity pro- 
file technique, while the cross marks the estimated origin of the 
event. White (black) is an increase (decrease) in intensity from 
the base image. 

The ratio of the BD image and the pre-event image was then cal- 
culated, giving a PBD image. This is described by the equation, 



h-k 



x 100, 



(1) 



where /, is the image at any time t and Iq is the pre-event im- 
age. This technique produces images that highlight the CBF and 
any associated dimming regions, with the intensity values of 
any given pixel corresponding to the percentage change in in- 
tensity with respect to th e pre-event image (for more details see 
IWills-Davev et alJl2007h . 

Figure Q] shows PBD images for the four events studied here 
from STEREO-B, with the event dates indicated in the lower 
right of each panel. The erupting AR was close to disk centre 
for both spacecraft for both the 2007 May 19 and 2007 Dec 07 
events due to their small separation. However the eruption was 
close to disk centre for STEREO-B but on the limb for STEREO- 
A for both the 2009 Feb 12 and 2009 Feb 13 events. As a result, 
only the STEREO-B observations are considered here. The ob- 
served pulse in each case extended over a large fraction of the 
solar disk and was observed in multiple images from both pass- 
bands. The observing cadence of the 195 A passband was 600 s 
for each event, while the 171 A passband operated at a cadence 
of 150 s for each event with the exception of the 2009 Feb 13 
event where it had a cadence of 300 s. 

3. Methods 

It is well known that the use of point-and-click techniques in 
conjunction with PBD images results in large errors that are 
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Fig. 2. PBD image intensity profiles 
(crosses) for the 2007 May 19 event as 
obtained from STEREO- A in the 195 A 
(left) and 171 A (right) passbands. The 
Gaussian fit (solid curve) to the positive 
section of each profile has been overplotted 
on each panel. The time of the leading image 
(/,) in each case is on the upper right of 
the panel. More cases are displayed in the 
online Figures [3]to[7] 



unique to each user. As a result, we developed an algorithm to 
automatically identify the presence of CBFs in PBD images and 
produce estimates of the kinematics of these CBFs. This algo- 
rithm is outlined in Section 13.11 with the resulting kinematics 
presented in Section l3~2l The algorithm was also used to investi- 
gate the variation in both the width (Section [3~3l and position an- 
gle (PA) -averaged integrated intensity (Section l3~4t of the CBF 
pulse. 



3.1. Pulse Identification 

The source of the CBF disturbance was identified from the PBD 
images by fitting ellipses to the visually traced wavefronts from 
the first two observations of the disturbance in both the 171 A 
and 195 A passbands, with the mean of the centres of the four 
ellipses taken as the origin of the CBF and the standard devi- 
ation giving the associated error (typically ~20 Mm). Next, a 
great-circle sector (i.e., an area on the sphere bounded by two 
great circles) projected onto the Sun was identified in which the 
disturbance was clearly visible throughout its propagation, with 
the source of the disturbance acting as the crossing point of the 
two bounding great circles (see Figure [1}. From each image the 
PBD intensity within this sector was then averaged across PA, in 
annuli of increasing radii with 1 degree width on the surface of 
the sphere, to produce an intensity profile as a function of dis- 
tance awa y from the source location (cf. similar techniques pro- 
posed by Warmuth et al.l 2004bl IPodladchikova & Berghmansl 
120051 1 Wills-Davevll2006l: I Veronig et al .1120101) . This process was 
repeated to create BD intensity profiles as well. 

The intensity profile produced by this technique was fit- 
ted using a Gaussian model because CBFs were noted by 



IWills-DavevI ( 12006b to have a Gaussian form. The technique 
may be viewed in the attached movie (see the movie attached 
to the online Figure |H), while the resulting intensity profiles for 
the 2007 May 19 event as obtained by STEREO- A in the 171 A 
and 195 A passbands are displayed in Figure [2] Note that the 
first three observations in the 195 A passband and the first six 
observations in the 171 A passband from both spacecraft were 
used for this work. Beyond this, the \ 2 values of the Gaussian 
fits were too large for accurate analysis. 

The intensity profile technique is judged to be more accurate 
than a normal point-and-click technique as it is semi-automated 
and reproducible, with any associated errors quantifiable as they 
result from the fitting of a Gaussian to an intensity profile across 
the pulse. In contrast, the point-and-click methodology com- 
monly used to identify CBFs is highly user-dependent, with 
large unquantifiable errors in the identification of the pulse posi- 
tion. Although a pulse may be tracked over larger distances using 
a visual method, increased image processing is often required to 
improve the visibility of the pulse. However, the intensity pro- 
file technique used here applies the same processing for all im- 
ages, with the result that measurements of the pulse are directly 
comparable. The intensity profile technique can also be used to 
process large amounts of data (such as those from SDO) rapidly 
as it does not require the same degree of user input as the visual 
method. 

In comparing different analysis techniques, it was also noted 
that BD and PBD images produced a more accurate estimate 
of the location of the pulse than running-difference images as 
they highlight relative motion of solar features with respect to a 
pre-event image rather than the previous image. This means that 
the bright pulse indicates the entirety of the CBF, rather than 
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Fig. 9. Distance-time plots for STEREO-A (top) and STEREO-B 
(bottom). The 171 A (circles) and 195 A (triangles) data have 
been combined as they follow similar kinematical curves. The 
mean offset distance, initial velocity, and acceleration terms re- 
sulting from the bootstrapping analysis (fit indicated by dashed 
line) are also stated in the bottom right of each panel, with errors 
represented by the standard deviation. The errors on each point 
are given by the error on the mean of the Gaussian fit applied 
to the intensity profile. More cases are displayed in the online 
FiguresHOlandQI] 

the portion of the pulse which has moved beyond its previous 
extent (for a discussio n of the prob lems associated with running- 
difference images see Attrill 2010). 

3.2. Pulse Kinematics 

The kinematics of the CBF pulse were studied by plotting the 
variation with time in the centroid of the Gaussian model applied 
to the positive section of the intensity profiles (see Figure|2). The 
errors associated with each data-point are given by the error as- 
sociated with the mean of the Gaussian fit to the intensity profile 
in each case. Data from both the 171 A and 195 A passbands 
were combined in this case as they have been observed to follow 
simila r kinematical curves dWarmufh et al.ll2004at IVeronig et al.l 
120101) . This was carried out for both STEREO spacecraft (see 
Figure [9} with the resulting plots fitted using a quadratic model 
of the form, 



r(f) = r + v t + -at, 



(2) 



where ro is the offset distance, vo is the initial velocity, a is the 
constant acceleration of the pulse, and t is the time elapsed since 
the first observation of the disturbance. The data were analysed 
using a residual resampling bootstrapping technique to ensure 
that the results were as statistically rigorous as possible. This 
technique works by fitting the given data (y,-; i = 1 , 2, . . . , n) with 
a specified model, yielding the fitted values ($) and residuals 
(e = yi - yi). The residuals are then randomly ordered, randomly 
assigned a sign (1 or -1), and added to the original fit values to 
produce a new data set. These data are fit using the same model 
with the resulting re-fitted parameters recorded. The process of 
randomizing residuals, applying them to the original fit and re- 
fitting is then repeated a large number of times (~ 10 000). This 
technique is statistically rigorous and produces a more accurate 



result than a simple model-fit to the given data (lEfronlll979l) . 
allowing each fit parameter to be characterised by a distribution. 

The inset values in Figure [9] are the mean and standard de- 
viation for these bootstrapped parameter distributions, which 
suggest that the pulse has an initial velocity of ~446 km s 
with a negative acceleration observed in both spacecraft (a 
-260 ± 130 m s~ 2 in STEREO-A and a -260 ± 150 m s~ 2 in 
STEREO-B). The errors associated with the acceleration terms 
are quite large (in each case the acceleration is negative within 
the 1-sigma error range). However, the kinematics are consistent 
between both STEREO spacecraft, suggesting that the decelera- 
tion is real. 

The initial velocity va lues given here ar e si milar to previ- 
ous es timates obtained by lLong et all (120081) and IVeronig et al.l 
(120081) who both studied this event, although in both of those 
cases a three-point Lagrangian interpolation technique was used 
to determine the kinematics of the CBF. Although that technique 
retains all of the data points through the use of an interpola- 
tion method, it has been observed to introduce artificial trends 
through the skewing of the interpolated edge points. As a result, 
it is possible to misinterpret the derived velocity and accelera- 
tion plots. In contrast, the bootstrapping technique used here is 
designed to determine the best-fit of a model to a given small 
data set, producing accurate estimates of the model parameters 
and quantifiable associated errors. 

3.3. Pulse Width 

The variation with dist ance of the full width at half maximum 
(FWHM; Ar = 2 V21n2cr) of the Gaussian fit applied to the PBD 
intensity profile was studied to identify any evidence of pulse 
broadening. The top two panels in Figure [12] show this variation 
for STEREO-A (left) and STEREO-B (right), with measurements 
from the differing passbands indicated by the different symbols. 
It is clear from these plots that an increase in pulse width with 
distance is present for both passbands as observed by both space- 
craft. This is indicative of pulse broadening and shows that the 
CBF spreads out spatially as it propagates. 

The variation in the temporal width of the pulse with time 
was also considered. The temporal width of the pulse, At, is 
defined as, 



Ar 



Ar 

Spills e 



(3) 



where Ar is the spatial width of the pulse and v pu i se is the ve- 
locity of the pulse. While an increase in the spatial width of 
the pulse may indicate dispersion, this could be negated by in- 
creasing pulse velocity, producing a pulse with constant tem- 
poral width. By examining the variation of the temporal width 
of the pulse with time it is possible to determine if the pulse 
is indeed dispersive. The resulting plot of the variation in pulse 
temporal width with time is given in the bottom two panels of 
Figure [L2l 

The temporal width of the pulse is also observed to increase 
with time, indicating that the pulse broadens in both space and 
time with propagation. The increase in spatial and temporal 
pulse width is apparent in the data from both STEREO space- 
craft, suggesting that it is a true feature of the disturbance and 
not an observational artifact. The rate of change of the spatial 
and temporal width of the pulse with distance and time respec- 
tively (d(Ar)/dr and d(AT)/dt) are given in TableQ] 

The broadening of the CBF pulse in b oth space and time con- 
firms the previous observations of both IWarmuth et al.l (1200 lb 
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Fig. 12. Top: Variation in pulse spatial width (Ar) with distance for the 2007 May 19 event. Bottom: Variation in pulse temporal 
width (At) with time. Panels contain data from STEREO-A (left) and STEREO-B (right). Pulse spatial width here refers to the 
FWHM of the fitted Gaussian pulse (i.e., Ar = 2 V21n2cr). The dashed line in all panels indicates the best linear fit to the data. More 
cases are displayed in the online Figures[T3lto[T5l 



and IVeronig et al.1 d2010h . These observations, when taken to- 
gether with the derived kinematics, are consistent with the inter- 
pretation of a CBF as a dispersive, decelerating pulse. 

3.4. Integrated Intensity 

The variation in the PA-averaged integrated intensity of the pulse 
was also examined in order to determine the physical nature of 
the disturbance. This was found by determining the variation in 
peak intensity with distance ((/Bo)max versus r) and the variation 
in the FWHM of the pulse with distance (Ar versus r). The PA- 
averaged integrated intensity could then be defined as, 

7 tot = Ar</ BD >max, (4) 

and plotted as a function of distance (see Figure [16] for the re- 
sulting plots for the 2007 May 19 event). 

This approach was used as the pulse has shown broadening 
(i.e., dispersion), with the result that the variation in peak inten- 
sity with distance may be influenced by this dispersion. The peak 
intensity of the BD intensity profile was used as this shows the 
actual emission of the pulse, while the PBD intensity shows the 
ratio of the pulse emission with respect to the background. 

Figure[l6]shows the variation in peak pulse intensity, FWHM 
and PA-averaged BD integrated intensity with distance for each 
spacecraft. The top plot shows a large discrepancy between the 
171 A and 195 A observations, with the 195 A data show- 
ing a large decrease with time while the 171 A data shows no 
strong variation. The FWHM does follow a linearly increasing 
trend with distance for both passbands, but the resulting varia- 
tion in PA-averaged BD integrated intensity is inconclusive in 
each case. In both cases the 195 A data appears flat with dis- 
tance (albeit with some stronger point-to-point variation), while 
the 171 A data shows a generally increasing trend. The rate of 
change of I m with distance (i.e., d(/ fof )/dr) is given in TableQ] 



These measurements show that more analysis is required to 
understand the morphology of the pulse. The higher cadence ob- 
servations available from SDO should allow the variation in peak 
intensity and PA-averaged integrated intensity with distance to 
be determined with a much higher degree of accuracy than pos- 
sible using STEREO. 



4. Results 

The intensity profile technique outlined in Section|2]is extremely 
effective at identifying the CBF pulse in observations, with the 
associated errors much lower than those achieved with previous 
techniques. Once the sector into which the pulse propagates has 
been identified, the technique automatically returns the centroid 
and FWHM of the fitted Gaussian as well as the integrated pulse 
intensity averaged over position angle at any given image time. 
The technique also produces accurate and reproducible estimates 
of the pulse characteristics, making these results more robust. 
The findings from the event of 2007 May 19 are discussed in 
detail in Section |4~T1 with several other CBF events summarized 
in Section l4~2l while all of the results are presented in TableQ] 

4.1. 2007 May 19 Event 

The event of 2007 May 19 displays kinematics consistent with 
those of a decelerating pulse, with an initial velocity that is to- 
wards the upper end of previous estimates. The acceleration term 
was observed to be negative within the 1-sigma error range, with 
the similar results from both spacecraft suggesting that these are 
the true kinematics of the pulse. 

The pulse was observed to display clear spatial and temporal 
broadening during propagation, indicating that the CBF is dis- 
persive. As a result, the variation with distance of the integrated 
intensity of the pulse averaged across position angle rather than 
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Fig. 16. Top: Variation in peak percentage base difference intensity ((/Bo)max) with distance. Middle: Variation in FWHM of the 
Gaussian pulse with distance. Bottom: Variation in integrated intensity (/,„,) with distance. Left-hand plots correspond to STEREO- 
A data, with STEREO-B data on the right-hand side. The dashed (dotted) line in all panels indicates the best linear fit to the 195 A 
(171 A) data. More cases are displayed in the online Figures [T7lto[T8l 



the peak intensity of the pulse was examined to try and under- 
stand the physical nature of the pulse. A decrease in the peak 
amplitude of a dispersive pulse does not necessarily imply a dis- 
sipative pulse (due to the presence of pulse broadening) and may 
result in a misinterpretation of the true nature of the disturbance. 

This event shows inconclusive variation in the PA-averaged 
integrated pulse intensity with distance from both observed pass- 
bands. The higher cadence 171 A data exhibits a generally in- 
creasing trend with distance, while the lower cadence 195 A 
observations show negligible variation on average, but strong 
point-to-point variation. The multi-passband, high-cadence ob- 
servations afforded by the Solar Dynamic Observatory (SDO) 
will allow the true variation (if any) to be determined to a high 
degree of accuracy. 

4.2. Further CBF Events 

Three additional CBF events from 2007 December 12, 
2009 February 12, and 2009 February 13 were also studied us- 
ing the intensity profile technique (see Figures |4]to |7j. The re- 
lationships discussed in Section [3] were plotted for these addi- 
tional events, with approximately similar results observed for 



each event. Table Q] shows the kinematics of all the pulses stud- 
ied and indicates that they displayed similar initial propagation 
velocities (~240-450 km s ). The 2007 May 19 event appears 
to have been a relatively fast event, with an initial velocity of 
~450 km s _1 and a statistically significant non-zero accelera- 
tion. The event of 2007 December 07 was much slower, with 
an initial velocity of ~260 km s and a statistically signifi- 
cant negative acceleration as observed by both STEREO space- 
craft. The quadrature events of 2009 February 12 and 13 were 
different from each other despite originating from the same re- 
gion. The 2009 February 12 event showed a faster initial veloc- 
ity (~405 km s ) and stronger deceleration (~ -291 m s~ 2 ), 
while the 2009 February 13 event had a slower initial veloc- 
ity (~274 km s -1 ) and a much weaker negative acceleration (~ 
-49 m s~ 2 ). The large errors associated with the acceleration 
terms given here indicate the difficulties associated with accu- 
rately determining the kinematics of CBFs from low cadence 
observations, despite the minimization of errors through the use 
of both the intensity profile technique and bootstrapping analy- 
sis. The distance-time plots for each event are given in the online 
Figures [10] and QT] for comparison. 
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Table 1. Wave properties of studied CBFs. 



Kinematics Dispersion Integrated Intensity 



Event 


Spacecraft 


Vo 




a 




d(Ar)/dr 


d(Ar)/dt 


dd,„ 


)ldr 






km s 


-l 


m s" 


2 






171 A 


195 A 


2007 May 19 


Ahead 


444 ± 


75 


-256 ± 


134 


0.3 ±0.1 


1.0 ±0.2 


1.08 ±0.53 


0.13 ± 1.07 




Behind 


447 ± 


87 


-260 ± 


149 


0.5 ±0.1 


1.1 ±0.2 


1.83 ±0.41 


0.24 ± 0.61 


2007 Dec 07 


Ahead 


270 ± 


37 


-152 i 


:73 


0.4 ±0.1 


1.0 ±0.3 


0.71 ± 0.37 


-0.049 ± 0.002 




Behind 


247 + 


62 


-117 ± 


110 


0.2 ± 0.2 


0.8 ± 0.3 


2.29 ± 0.41 


-1.07 ±0.12 


2009 Feb 12 


Behind 


405 ± 


93 


-291 ± 


166 


0.4 ± 0.2 


2.0 ± 0.5 


-3.27 ± 2.39 


-0.59 ±0.11 


2009 Feb 13 


Behind 


274 ± 


53 


-49 ± 


34 


0.4 ±0.1 


0.5 ±0.1 


4.33 ± 0.66 


-0.74 ± 0.36 



Notes. Kinematic values refer to the mean and standard deviation error of the bootstrapping parameter distributions. Dispersion and integrated 
intensity values refer to the rate of change of the relevant parameters, resulting from the linear fits shown in Fieures [T2land[76l 



Table[TJalso shows the rate of change of spatial and temporal 
pulse width Ar and At with distance and time respectively, as 
well as the rate of change of the PA-averaged integrated pulse 
intensity, AI, or , associated with each event (i.e., the slope of the 
lines shown in FiguresfT2land[T6Ti. The results indicate that all of 
the observed events exhibited clear pulse broadening in both the 
spatial and temporal domains, with the d(Ar)/dr and d(AT)/dt 
parameters positive for both the 171 A and 195 A passbands. 
This suggests that pulse broadening is a general characteristic of 
CBFs and must be accounted for by all theories that seek to ex- 
plain this phenomenon. The plots showing the variation in both 
spatial and temporal pulse width with distance and time respec- 
tively for each event are given in the online Figures [T3"lto [TBI 

The variation in PA-averaged integrated pulse intensity 
with distance is more difficult to interpret. Both of the online 
Figures [17] and [18] show the variation in peak intensity (top), 
FWHM (middle) and PA-averaged integrated pulse intensity 
(bottom) with distance. In each case, the 195 A peak intensity 
drops with distance while the 171 A peak intensity variation is 
inconclusive, although the FWHM of each passband tends to 
increase. For the 2007 December 07 event, the resulting PA- 
averaged integrated pulse intensity shows a slight decrease with 
distance for the 195 A data, while showing an apparent increase 
with distance for the 171 A data, with a large separation be- 
tween the 171 A and 195 A passbands. In contrast, both the 
171 A and 195 A observations tend to drop with distance for the 
2009 February 12 event, while the 2009 February 13 event shows 
an increase and decrease with distance for the 171 A and 195 A 
data respectively. These observations suggest that the variation 
in PA-averaged integrated pulse intensity with distance is not 
well-defined and requires further investigation. 



5. Conclusions 

Several CBFs were studied by applying a semi-automated 
technique to data from the EUVI telescopes onboard both 
STEREO spacecraft. This technique applies Gaussian fits to in- 
tensity profiles across propagating CBFs, returning the posi- 
tion, width, and PA-averaged integrated intensity of the pulse 
for each observation. This intensity profile approach is sim- 
ilar to CBF id e ntificat io n techniques previously prop o sed b y 



Warmutheta.1 (l2004bh: iPodladchikova & Berghmansl (120051) : 
Wills-Davevl ( 20061) : Veronig et al. d2010l) . However, by combin- 
ing it with a statistically rigorous bootstrapping method and the 
high cadence observations afforded by STEREO/EWl we can 
minimise the errors typically encountered with the analysis of 
CBFs. 



The similarity in derived velocity between this work and 
previous investigations is interesting given that most previous 
works have used point-and-click techniques applied to running- 
difference images. These studies identify the forward edge of the 
CBF at a given time, which is then used to determine the kine- 
matics of the disturbance as a whole. The analyses performed 
using such techniques have mainly returned kinematics that sug- 
gest a zero acceleration (i.e., constant velocity) interpretation 
of the CBF phenomenon. In contrast, our semi-automated tech- 
nique uses the pulse centroid derived from percentage base dif- 
ference intensity profiles to return a constant, non-zero acceler- 
ation (i.e., variable velocity) interpretation. In the presence of 
pulse broadening, as found here, the forward edge of a deceler- 
ating pulse will appear to move faster than the pulse centroid. 
This suggests that the true kinematics of CBFs may have previ- 
ously been disguised through use of the wrong position within 
the pulse profile to characterise its location. We also note that a 
variable acceleration may be present, but this can only be studied 
with an increased number of data points. 

A positive increase in the variation of the width of the 
Gaussian fit to the CBF pulse in both the spatial and temporal 
domains was observed for all events studied. While the increase 
in spatial pulse width is suggestive of a dispersive pulse, this is 
confirmed by the observed increase in the temporal width of the 
pulse. From Equation [3] an increase in the spatial pulse width 
may be negated by an increase in pulse velocity, producing a 
pulse with a constant temporal width. However, Figures [12] to 
[TBI show that this is not the case. This indicates that all of the 
CBFs studied showed significant pulse broadening and may be 
interpreted as being dispersive. Although some previous studies 
have suggested a dispersive nature for "EIT waves", the true ex- 
tent of this dispersion may have been disguised by the nature of 
point-and-click analyses of running-difference images. 

The variation in PA-averaged integrated pulse intensity with 
distance was studied rather than the variation in peak pulse in- 
tensity due to the presence of pulse broadening. This was found 
by multiplying the FWHM by the peak pulse intensity at a given 
time, and produced inconclusive results. While the peak inten- 
sity of the pulse was generally observed to decrease with dis- 
tance and the FWHM was observed to increase with distance 
in both passbands, the PA-averaged integrated intensity typi- 
cally showed no strong variation with distance in either pass- 
band, although a strong offset between the 171 A and 195 A 
data was typically observed. Contrasting results were found for 
the 2009 February 12 and 2009 February 13 events respectively, 
suggesting that further investigation using simultaneous analysis 
of multiple passbands at very high cadence is required, some- 
thing that will be routinely available from the SDO spacecraft. 
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The results of our analysis suggest that the studied CBFs may 
be interpreted as dispersive pulses exhibiting negative accelera- 
tion. This is consistent with the fast-mode magnetoacoustic wave 
interpretation for a freely-propagating "EIT wave". The varia- 
tion in the integrated intensity of the pulse was inconclusive, im- 
plying that more analysis is required to definitively determine 
the physical nature of the disturbance. 

The consistency of the results between the events studied 
suggests that the conclusions drawn here may be applicable to 
a larger sample of CBFs. The initial velocities of the CBFs 
are compara ble to the lower range of estimated Alfven speeds 
proposed by IWills-Davev et al.l (120071) . suggesting that the ran- 
domly structured nature of the coronal magnetic field may have 
an important effect on the pro pagation of CBFs; an e ffect shown 
in simulations performed by iMurawski et al.l (1200 lb . Although 
the kinematics and dispersive nature of the pulses suggest a 
magnetoacoustic wave interpretation may be correct, further 
study is required. The ~10 s temporal cadence across multiple 
EUV passbands available with the launch of the SDO spacecraft 
should allow these results to be studied in much greater detail. 
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Fig. 3. Same as for Figure [2] but for 
2007 May 19 STEREO-B. 
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Fig. 5. Same as for Figure [2] but for 
2007 Dec 07 STEREO-B. 
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Fig. 7. Same as for Figure [2] but for 
2009 Feb 13 STEREO-B. 
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Fig. 8. Left panel: PBD image with arc sector over-plotted. Right panel: Intensity profile across arc sector. The intensity profile at a 
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Fig. 10. Same as Figure|9]but for the event on 2007 Dec 07. 
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13. Same as for Figure[T2]but showing the event on 2007 Dec 07. 
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14. Similar to FigureQ/J] now only showing STEREO-B data for the event on 2009 Feb 12. 



Long et al.: Wave properties of CBFs, Online Material p 6 

STEREO-B 





300 




250 


? 


200 


s 






150 


In 




< 


100 




50 








A 195 A 
O 171 A 



A 



O 



O 



.0.. 



1500 

1000 

500 




100 200 300 400 500 600 
Distance (Mm) 

A ' 1 



A 


195 


A 





171 


A 



A 



1-® ffl 

D ..................... , 

05:40 05:50 06:00 06:10 06:20 
Start Time (13-Feb-09 05:35:31) 



Fig. 15. Similar to Figure Q~2] now only showing STEREO-B data for the event on 2009 Feb 13. 
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Fig. 17. Same as for Figure[T6]but for the event of 2007 December 07. 
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Fig. 18. Similar to Figure Q2] but showing the STEREO-B data from the events of 2009 February 12 (left) and 2009 February 13 
(right). 



